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UNBOUNDED SUPERSOLUTIONS OF SOME 
QUASILINEAR PARABOLIC EQUATIONS: A 
DICHOTOMY 

JUHA KINNUNEN AND PETER LINDQVIST 


Abstract. We study unbounded ’’supersolutions” of the Evolu¬ 
tionary p-Laplace equation with slow diffusion. They are the same 
functions as the viscosity supersolutions. A fascinating dichotomy 
prevails: either they are locally summable to the power p—1 + ^ — 0 
or not summable to the power p — 2. There is a void gap between 
these exponents. Those summable to the power p — 2 induce a 
Radon measure, while those of the other kind do not. We also 
sketch similar results for the Porous Medium Equation. 


1. Introduction 

The unbounded supersolutions of the Evolutionary p-Laplace Equa¬ 
tion 

I 

— - V-(|Vu|P-Vu) =0, 2<p<oo, 

exhibit a fascinating dichotomy in the slow diffusion case p > 2. This 
phenomenon was discovered and investigated in jTT]. The purpose of 
the present work is to give an alternative proof, directly based on the 
iterative procedure in [Sj. Besides the achieved simplification, our proof 
can readily be extended to more general quasilinear equations of the 
form 

du „ ^ „ 

— - V ■ A(a;, f, M, Vn) = 0, 

which are treated in the book [DGV]. The expedient analytic tool is 
the intrinsic Harnack inequality for positive solutions, see [6]. We can 
avoid to evoke it for siipersolutions. We also mention the books |1] and 
m as general references. 

The supersolutions that we consider are called p-supercaloric func- 
tion£|. They are pointwise defined lower semicontinuous functions, fi¬ 
nite in a dense subset, and are required to satisfy the Comparison 
Principle with respect to the solutions of the equation; see Definition 
12.41 below. The definition is the same as the one in classical potential 
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^They are also called p-parabolic functions, as in [S]. 
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theory for the Heat EquatioiJl, to which the equation reduces when 
p = 2, see 12^. Incidentally, the p-supercaloric functions are exactly 
the viscosity supersolutions of the equation, see 

There are two disjoint classes of p-supercaloric functions, called class 
IB and 971. We begin with the former one. Throughout the paper we 
assume that is an open subset of M” and we denote Qt = x (0, T) 
for T > 0. 


Theorem 1.1 (Class IB). Let p > 2. For a p-supercaloric function 
V : Qt (—cxo, oo] the following conditions are equivalent: 

(i) ^ e L\;^{ytT), 

(ii) the Sobolev gradient Vn exists and Vn G whenever 

q' + 

(hi) V e whenever q < p — 1 + ^. 


In this case there exists a non-negative Radon measure p such that 

( 1 - 2 ) r 


'0 JO . 


-I-(|Vn|^ ^Vv,Vip}'^ dxdf = J p dp 


for all test functions p G ( 7 ^( 127 ’). In other words, the equation 

‘L - V ■ (\vvr^vvp 


holds in the sense of distributions, cf. na. It is of utmost importance 
that the local summability exponent for the gradient in (ii) is at least 
p — 2. Such measure data equations have been much studied and we 
only refer to [2]. For potential estimates we refer to ng, ng. 

As an example of a function belonging to class 93 we mention the 
celebrated Barenblatt solution 


(1,3) 



r n 


93(x,f) = < 




. 0, 

when f < 0, 


when f > 0, 


found in 1951, cf. |T]. Here A = n{p — 2) +p and p > 2. It is a solution 
of the Evolutionary p-Laplace Equation, except at the origin x = 0, 
t = 0. Moreover, it is a p-supercaloric function in the whole R"' x R, 
where it satishes the equation 

(993 

— - V - (|V93|^’-2V93) = cS 

in the sense of distributions (5 = Dirac’s delta). It also shows that the 
exponents in (i) and (ii) of the previous theorem are sharp. 


^Yet, the dichotomy we focus our attention on, is impossible for the Heat 
Equation. 
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A very different example is the stationary function 

j 

where the g^’s are an enumeration of the rationals and the Cj > 0 are 
convergence factors. Indeed, this is a p-supercaloric function, it has a 
Sobolev gradient, and v{qj,t) = oo along every rational line x = qj, 
—oo < t < oo, see [T8] . 

Then we describe class 


Theorem 1.4 (Class VJt). Let p > 2. For a p-supercaloric function 
V : Qt oo] the following conditions are equivalent: 

(i) v<tpz\aT), 

(ii) there is a time to, <to <T, such that 


liminf v{y,f){t — tQ)p-^ > for all x e Ll. 

t>to 


Notice that the inhnities occupy the whole space at some instant to- 
As an example of a function from class DJI we mention 

f 1 

-i—, when t > to, 

(t-to)^ 

I 0, when t < to, 

where il G (7(11) fl lT)}’^(f2) is a weak solution to the elliptic equation 

V ■ (|Vit|^-'Vit) +^it=0 

and il > 0 in fl. The function QJ is p-supercaloric in fl x M, see equation 
fl3.3p below. This function can serve as a minorant for all functions 
T > 0 in 911. No (T-hnite measure is induced in this case. As far as we 
know, these functions have not yet been carefully studied. 

A function of class DJI always affects the boundary values. Indeed, at 
some point on (^O) to) on the lateral boundary dQ x (0, T) it is necessary 
to have 

limsup v{x,t) = oo. 

(x,t)->((o,to) 

This alone does not yet prove that v would belong to 91t. A convenient 
sufficient condition for membership in class 93 emerges: If 


limsup v{x,t) < oo for every (^,t) G dfl x (0,T), 

(x,T)^(^,t) 

then V ^ DB. 

It is no surprise that a parallel theory holds for the celebrated Porous 
Medium Equation 

dll 

— -A('u™') = 0, 1 < m < cxo. 

dt ^ ^ ’ 
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We refer to the monograph |22] about this much studied equation. We 
sketch the argument in the last section. 


2. Preliminaries 


We begin with some standard notation. We consider an open domain 
12 in R”' and denote by LP{ti,t 2 ] the Sobolev space of functions 

V = v{x,t) such that for almost every t,ti < t < t 2 , the function 
X H-)■ v{x,t) belongs to and 



where Vu = (^, • • • , -^) is the spatial Sobolev gradient. The dehni- 

tions of the local spaces U’{ti,t 2 ] ( 12 )) and t 2 ] are 

analogous. We denote = 12 x ( 2 i,t 2 ) and recall that the parabolic 
boundary of khe set (12 x {ti}) U (<912 x ( 2 i, 22 ))- 

Definition 2.1. A function u E L^(ti,t 2 ; IV^’^(Q)) is a weak solution 
of the Evolutionary p-Laplace Equation in 124^ ,* 2 ; ik 



( 2 . 2 ) 


for every 99 G C^{Qti,t 2 )- If) in addition, u is continuous, then it 
is called a p-caloric function. Further, we say that u is a weak su¬ 
persolution, if the above integral is non-negative for all non-negative 
(p E C^{Qti,t 2 )- If fh® integral is non-positive instead, we say that u is 
a weak subsolution. 

By parabolic regularity theory, a weak solution is locally Holder 
continuous after a possible redehnition in a set of n -|- 1 -dimensional 
Lebesgue measure zero, see ESI and jl]. In addition, a weak superso¬ 
lution is upper semicontinuous with the same interpretation, cf. ca. 

Lemma 2.3 (Comparison Principle). Assume that 


u, V E L^[ti,t 2 ; W^’^(12)) n (7(12 x [^ 1 , 22 ))- 


If V is a weak supersolution and u a weak subsolution in Ilti,t 2 such 
that V > u on the parabolic boundary ofQti,t 2 ! then v > u in the whole 



The Comparison Principle is used to dehne the class of p-supercaloric 
functions. 

Definition 2.4. A functions : 12 t^^i 2 (—cx 3 ,cx 3 ] is called p-superca/onc. 


if 


(i) V is lower semicontinuous, 

(ii) V is hnite in a dense subset. 
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(iii) V satisfies the comparison principle on each interior cylinder 

(E If ^ is a p-parabolic function in 

and if h < u on the parabolic boundary of then 

h < u in the whole Df/ +/. 

— ^ 1,^2 

We recall a fundamental result for hounded functions, which is also 
applicable to more general equations. 

Theorem 2.5. Let p > 2. If v is a p-supercaloric function that is 
locally hounded from above in Qt, then the Soholev gradient Vv exists 
and Vv G Moreover, v G T; Wlff{VL)) and v is a weak 

supersolution. 

A proof based on auxiliary obstacle problems was given in |2], The¬ 
orem 1.4. A more direct proof with infimal convolutions can be found 
in [T7] . 

In order to apply the previous theorem, we need bounded functions. 
The truncations 


Vj{x,t) = mm{v{x,t),j}, j = 1,2,..., 

are p-supercaloric, if v is, and since they are bounded from above, they 
are also weak supersolutions. Thus Vvj is at our disposal and estimates 
derived from the inequality 

where > 0 and ip G sxe available. The starting point for 

our proof is the following theorem for the truncated functions. 

Theorem 2.7. Let p > 2. Suppose that v > 0 is a p-supercaloric func¬ 
tion in Qt with initial values v{x, 0) = 0 m Q. Ifvj G Lp( 0, T; ITQ^’^(r2)) 
for every j = 1,2,..., then 

(i) V G L'^{flTi) whenever q < p — 1 ^ and Ti < T, 

(ii) the Sobolev gradient Vu exists and Vu G L'^'(r2ri) whenever 
q' < p — 1 + and Ti < T. 

Proof. See [^. □ 

We remark that the summability exponents are sharp. It is decisive 
that the boundary values are zero. The functions of class HOI cannot 
satisfy this requirement. As we shall see, those of class HB can be 
modified so that the theorem above applies. 

The standard Caccioppoli estimates are valid. We recall the following 
simple version, which will suffice for us. 
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Lemma 2.8 (Caccioppoli). Let p >2. If u > 0 is a weak subsolution 
in Qt, then the estimate 



+ esssup / f{xyu{x,tY dx 
2 ti<t<t2 Jn 

rt2 r r 


<C{p) 


uP\VC\Pdxdt+ / Cixfuix^ty 


Iti Jq 


fn 


^2 


da; 




holds for every ( = ({x) > 0 in 0 < ti < t 2 < T. 


Proof. A formal calculation with the test function cj) = vCy gives the 
inequality. See [1], |9]. □ 


Infimal Convolutions. The inhmal convolutions preserve the p-supercaloric 
functions and are Lipschitz continuous. Thus they are convenient ap¬ 
proximations. If n > 0 is lower semicontinuous and hnite in a dense 
subset of VLt, then the infimal convolution 

vyx,t)= inf \v{y,T) + ^{\x-yy + \t-Ty)\ 

(y,T)£nT t ‘is J 

is well dehned. It has the properties 

• v^{x, t) v{x, f) as e —)■ 0, 

• is locally Lipschitz continuous in Qt, 

• the Sobolev derivatives ^ and exist and belong to 

Assume now that n is a p-supercaloric function in Qt- Given a subdo¬ 
main D (E Qt, the above v'^ is a p-supercaloric function in D, provided 
that e is small enough, see [H]. 


3. A Separable Minorant 

We begin with observations, which will simplify some arguments 
later. 

Extension to the past. If n is a non-negative p-supercaloric function 
in IIt’, then the extended function 

, , \v(x,t), when 0 < f < T, 

vix, t) = < 

I 0, when f < 0, 

is p-supercaloric in x (—oo,T). We use the same notation for the 
extended function. 

A separable minorant. Separation of variables suggests that there 
are p-caloric functions of the type 

v{x, t) = {t — tQ)~^u{x). 
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Indeed, if is a domain of finite measure, there exists a p-caloric 
function of the form 

(3.1) 5J(a:,t) =-when t > to, 

{t - to)t^ 

where it G C{^1) fl WQ^’^(r2) is a weak solution to the elliptic equation 

(3.2) V ■ + ^it = 0 

and il > 0 in n. The solution if is uniqu^. (Actually, il G for 

some exponent a = a{n,p) > 0.) The extended function 


(3.3) 


Tr(a;, t) 


if(x) 

{t - to)^ 
0 , when 


when 
t < to- 


t > to, 


is p-supercaloric in hi x M. The existence of il follows by the direct 
method in the Calculus of Variations, when the quotient 

(/^w^da:)" 


is minimized among all functions w in with w ^ 0. Replacing 

w by its absolute value |tc|, we may assume that all functions are non¬ 
negative. Sobolev’s and Holder’s inequalities imply 

J{w) > c(p, n)|r2|^“"“2, 


for some c{p,n) > 0 and so Jq = inf^ J{'w) > 0. Choose a minimizing 
sequence of admissible normalized functions Wj with 

lim J{wj) = Jo and ||u',||Lp(n) = 1- 

j^oo 

By compactness, we may extract a subsequence such that ^ Vw 
weakly in L^{VL) and Wj^ —)■ w strongly in LP{fl) for some function w. 
The weak lower semicontinuity of the integral implies that 


J{w) < lim inf Jiwj^) = Jo- 

k^oo 

Since w G this means that tc is a minimizer. We have tc > 0, 

and tc ^ 0 because of the normalization. 

It follows that w has to be a weak solution of the Euler-Lagrange 
equation 

V ■ + JoWw^-^^^w = 0 

with ||tc||Lp(r 2 ) = 1- By elliptic regularity theory w G ^(n), see [20] . 
Finally, since V ■ (|< 0 in the weak sense and tc > 0 


^Unfortunately, the otherwise reliable paper [J. GarcTa Azorero, I. Peral 
Alonso: Existence and nonuniqueness for the p-Laplacian: Nonlinear eigenvalues, 
Communications in Partial Differential Equations 12, 1987, pp. 1389-1430], con¬ 
tains a misprint exactly for those parameter values that would yield this fuirction. 
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we have that z/; > 0 by the Harnack inequality [20]. A normalization 
remains to be done. The function 

il = Cu, where JoC^~^ = 

will do. 


One dimensional case. In one dimension the equation is 

4 - (= 0, 0 < X < L. 

ax V / p — 2 

It has the hrst integral 


p — 1 
p 




il2 

2{p-l) 


C 


in the interval [0, L], Now if(0) = 0 = iI(T) and il'(f) = 0. This deter¬ 
mines the constant of integration in terms of il'(O) or of the maximal 
value M = maxil = Solving for if', separating the variables, and 

integrating from 0 to -I, one easily obtains the parameters 

M = Ci{p)L^ and it^(O) = —it'(L) = C 2 {p)L^. 


The constants can be evaluated. In passing, we mention that has 
interesting properties as a special function. 


4 . Harnack’s Convergence Theorem 


A known phenomenon for an increasing sequence of non-negative p- 
caloric functions is described in this section. The analytic tool is an 
intrinsic version of Harnack’s inequality, see [?], pp. 157-158, |5], and 


Lemma 4.1 (Harnack’s inequality). Let p > 2. There are constants 
C and 7 , depending only on n and p, such that if u > 0 is a lower 
seniicontinuos weak solution in 

H(xo, 4 i 7 ) X (to — 40 , to + 40 ), 
then the inequality 


where 0 = 




u{xo,to)P 2 ’ 


( 4 . 2 ) 


n(xo,to)<7 inf M(x,to + 0) 


is valid. 


Notice that the waiting time 0 depends on the solution itself. 

Proposition 4.3. Suppose that we have an increasing sequence 0 < 
hi < h 2 < h^ <... of p-caloric functions in Qt denote h = 
limfc^oo hfc- If there is a sequence {xk, tu) — t (xq, to) such that hk{xk, tk) 
-l-cxo, where xq G 12 and 0 < to < T, then 

liminf h{y,t){t — to)^ > 0 for all x G 12. 
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Thus, at time t^, 


lim h{y, t) = oo in 

t>to 

Remark 4.4. The limit function h may be finite at every point, though 
locally unbounded. Keep the function in mind. — The proof will 
give 


h(x, t) > 

{t - to)^ 

Proof: Let B{xq,4:R) <1 hi. Since 

CRP 


Ok — 


in n X {to, T). 


—>■ 0 , 


hk{xk,tk)P ^ 

Harnack’s Inequality fl4.2l) implies 

(4.5) hk{xk, tk) < ^hk{x, tk + 9k) 

when X G B{xk,R) provided B{xk,‘^R) x (4 —40^,4+ 46*^) d VLt- The 
center is moving, but since Xk —t xq, equation fl4.5p holds for sufficiently 
large indices. Let A > 1. We want to compare the solutions 

il^(x) 1 , / X 

- — and tik{x,t) 

(f - 4 + (A - 1)4) 

when t = tk + 9k and x G B{xo, R)- Here is the positive solution of 
the elliptic equation fl3.2l) in B{xo,R) with boundary values zero. We 
get 


il^(x) 


{t — tk + {R — 1)4) 

il^(x 


p-2 


il^(x) 




{KCR 


:p p-2 


1 hk {Xk , tk) 


< 


{ACRp)p-^ 

by taking A so large that 


—74(x, 4 + 4) < hk{x, tk + 9k) 


7||h’'II 


L°°(B{xo,R)) 


< 1. 

{ACRP)^ 

By the Comparison Principle 

- ^ - — < hk{x,f) < h{x,t) 

(t - 4 + (A - 1)4) 

when t >tk + 9k and x G B{xo, R). By letting k ^ oo, we arrive at 
il^(x) 


{t - to) 


< h{x, t) when to < t < T. 


p-2 
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Here it^ depended on the ball B{xq, R), bnt now we have many more 
inhnities, so that we may repeat the procednre in a snitable chain of 
balls to extend the estimate to the whole domain □ 


Proposition 4.6. Suppose that we have an increasing sequence 0 < 
hi < h 2 < /13 < ... of p-caloric functions in Qt denote h = 
limfc^oo hfc. If the sequence {hk\ is locally bounded, then the limit func¬ 
tion h is p-caloric in Qt- 

Proof. In a strict subdomain we have the Holder continuity estimate 

\hk{,xi,ti) - hk{x2,t2)\ < C\\hk\\ (^|X2 - Xi|“ + 1^2 - 

so that the family is locally equicontinuous. Hence the convergence 
hfe —)■ h is locally uniform in Qt. Theorem 24 in [LM] implies that 
{V/ifc} is a Cauchy sequence in Thus we can pass to the 

limit under the integral sign in the equation 


-hk^ + (iVhfcl^’-^Vhfc, V(p)^ dxdt = 0 


Jo Jn 

as k ^ 00 . From the Caccioppoli estimate 

[ C|V/ifc|Pdxdt 

Jti Jn 

<C{p) [ f hl\VC\^dxdt C{p) [ C{xYhk{x,tf 

Jti Jn Jn 

we deduce that h G Lf^^( 0 ,T; Wl^fiVt)). 


t2 


dx 


□ 


5. Proof of the Theorem 

For the proof we start with a non-negative p-supercaloric function 

V dehned in VLt- By the device in the beginning of Section [21 we £x 
a small h > 0 and redehne v so that v{x,t) = 0 when t < 5. This 
function is p-supercaloric. This does not affect the statement of the 
theorem. The initial condition x(x,0) = 0 required in Theorem 12.71 is 
now in order. 

Let Q 21 CC 12 be a cube with side length 41 and consider the con¬ 
centric cube 

Qi = {x||xi - x°| < l,i = 1,2,.. .n} 

of side length 21. The center is at x'^. The main difficulty is that 

V is not zero on the lateral boundary, neither does vj obey Theorem 
ED We aim at correcting v outside Qi x (0, T) so that also the new 
function is p-supercaloric and, in addition, satishes the requirements of 
zero boundary values in Theorem 12.71 Thus we study the function 

V in Qix{0,T), 
h in {Q 2 i\Qi) X {fl,T), 


(5.1) 


w = 
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where the function h is, in the outer region, the weak solution to the 
boundary value problem 

fh = 0 on dQ 2 i x (0,T), 

(5.2) i h = n on dQi x (0, T), 

[/i = 0 on {Q 2 i\Qi) X {0}- 

An essential observation is that the solution h does not always exist. 
This counts for the dichotomy. If it exists, the truncations Wj satisfy 
the assumptions in Theorem 12.71 as we shall see. 

For the construction we use the inhmal convolutions 

v%x,t)= inf \v{y) + ^{\x-y\‘^ + \t-T\‘^)\. 

{y,T)GQT t 2£ J 


They are Lipschitz continuous in Q 21 x [0,T] and weak supersolutions 
when e is small enough. Then we dehne the solution as in formula 
fl5.2p above, but with in place of v. Then we dehne 


w 


€ 


in Qi X (0,T), 
in {Q 2 i\Qi) X {0,T), 


and w^{x, 0 ) = 0 in Now < v^, and when t < S we have 0 < < 

= 0 so that h^{x,t) = 0 when t < 6. The function satishes the 
comparison principle and is therefore a p-supercaloric function. Here 
it is essential that < v^. The function is also (locally) bounded; 
thus we have arrived at the conclusion that is a weak supersolution 
in Q 21 X (0,T). 

There are two possibilities, depending on whether the sequence {h'^} 
is bounded or not, when £ \ 0 through a sequence of values. 

Bounded case. Assume that there does not exist any sequence of 
points (xgjfg) —)■ {xo,to) such that 


lim h^(x^, ts) = 00, 
£->•0 


where Xq G Q 21 \ Qi and 0 < to < T (that is an interior limit point). 
By Proposition 14.61 the limit function h = limg^o is p-caloric in its 
domain. The function w = limg^o itself is p-supercaloric and agrees 
with formula fl5.ip . 

By Theorem 12.51 the truncated functions Wj = min{i/;(x,t), j}, j = 
1,2,..., are weak supersolutions in Q 21 x (0, T). We claim that 

Wj e LP{ 0 ,r;W^'P{Q 2 i)) when T'< T. 

This requires an estimation, where we use 

L = snp{h{x, t) : {x, t) G {Q 21 \ Q 51 / 4 ) x (0, T')}. 

Let C = C(a;) be a smooth cutoff function such that 0<C<1)C = 1 
in Q 21 \ Q 31/2 and C = 0 in < 55 ^/ 4 . Using the test function when 
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deriving the Caccioppoli estimate we get 


'0 J Q2l\Q3l/2 


I'Vwjf dx dt 


< 


'0 J Q2l\Q3l/2 
rT' 


\Vh\^ dxdt < 


'O J Q2i\Qu/A 


C^|Vh|?’dxdt 


<C{p) 


'0 J Q2i\Qi 


hP\VC\Pdxdt + 


'Q2l\Q5l/4, 


h{x, T'Y dx 


< C{n,p){U>r-^T + 


where we used the fact that |Vtaj| = |Vmin{h,j}| < |Vh| in the 
outer region. Thus we have an estimate over the outer region Q21 \ 
Qzil2- Concerning the inner region Q31/2, we hrst choose a smooth 
cutoff function p = p{x,t) such that 0<i7<l, i7 = lin Q31/2 and 
r; = 0 in Q21 \ Qgi/i- Then the Caccioppoli estimate for the truncated 
functions Wj, j = 1,2,..., takes the form 


r-T' 


'0 •JQ 31/2 


<Cf 


dxdt < 


'0 J Q 21 


p'^\Vwj\^ dx dt 


\Vp\Pdxdt + Cf 


\ptY dx dt. 


'0 •JQ 21 


Jo JQ 21 

Thus we have obtained the estimate 

f [ \Vwj\Pdxdt<Cf 

Jo J Q 21 

over the whole domain Q21 x ( 0 , T') and it follows that Wj G L( 0 , T'; WY^{Q2i))- 
In particular, the crucial estimate 

/ / iVtcil^dxdf < 00, 

Jo J Q 21 

which was taken for granted in [ 10 ], is now established^ 

From Theorem 12.71 we conclude that v G L^{Qi) and Vn G {Qi) 
with the correct summability exponents. Either we can proceed like 
this for all interor cubes, or the following case occurs. 

Unbounded case. If there is a sequence (xo,to) such that 

lim h^ixe, tY) = 00 

e^O 

for some Xq G Q21 \ Qi, 0 < to < T, then 

v{x,t) > h{x,t) > (t — to)~^ii{x), 


^The class 911 passed unnoticed in m- 
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when t > to, according to Proposition 14.31 Thus v{x,to+) = oo in 
Q 21 \ Qi- But in this construction we can replace the outer cube with 
that is, a new h is dehned in \ Qz- Then by comparison 

v>h^> 

and so v{x, to+) = cxo in the whole boundary zone \ Qi. 

It remains to include the inner cube Qi in the argument. This is 
easy. Reflect h = in the plane xi = x\ +1, which contains one side 
of the small cube by setting 

h*{xi,X 2 - ,Xn) = h{2x\ + 21 - Xi,X 2 , . . . ,Xn), 

SO that 

xi + (2(a;? + /) - xi) 

-^^ = x\ + l 

as it should. Recall that x° was the center of the cube. (The same can 
be done earlier for all the k^.) The reflected function h* is p-caloric. 
Clearly, v>h*hy comparison. This forces v{x,to+) = 0 when x G Qi, 
Xi > x\. A similar reflexion in the plane xi = x\ — I includes the other 
half Xi < x^. We have achieved that v{x,to+) = 00 also in the inner 
cube Qi- This proves that 

v{x,to+) = 00 in the whole 


6. The Porous Medium Equation 


We consider the Porous Medium Equation 

du 


dt 


- A(n”^) = 0 


in the slow diffusion case m > 1. The equation is reated in detail 
in the book j22]- We also mention |21] and [11]. In CD the so- callecH 
viscosity supersolutions of the Porous Medium Equation were defined in 
an analoguous way as the p-supercaloric functions. Thus they are lower 
semicontinuous functions v : Qt [0,C)o], finite in a dense subset, 
obeying the Comparison Principle with respect to the solutions of the 
equation. 

Again we get two totally distinct classes of solutions, called class 
and DJI. Now the discriminating summability exponent is m — 1. We 
begin with 03. 


Theorem 6.1 (Class ®). Let m > 1 . For a viscosity supersolution 
V : Qt [0, cxo] the following conditions are equivalent: 

(i) n G LZ;\nT), 

(ii) the Sobolev gradient V(n™'“^) exists and G 

whenever q' < 1 + 

^The label ’Viscosity” was dubbed in order to distinguish them and has little to 
do with viscosity. The name ’’m-superporous function” would perhaps do instead? 
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(iii) V G whenever q < m + ^. 

A typical member of this class is the Barenblatt solution for the 
Porous Medium Equation. In this case a viscosity supersolution is a 
solution to a corresponding measure data problem with a Radon mea¬ 
sure in a similar fashion as for the Evolutionary p-Laplace Equation. 
The other class of viscosity supersolutions is 971. Unfortunately, this 
class was overlooked in m- 

Theorem 6.2 (Class 9Jt). Let m > 1. For a viscosity supersolution 
V : Qt [0, cxo] the following conditions are equivalent: 

(ii) there is a time to, <to <T, such that 

liminf v{y,t){t — to)^^^ >0 for all x G U. 

t>to 

Notice that again the inhnities occupy the whole space at some in¬ 
stant to- In m Theorem 3.2 it was established that a bounded vis¬ 
cosity supersolution n is a weak supersolution to the equation: n™' G 

whenever tp G (7“(fir) and if > 0. 

We shall deduce the above theorems from the following result. 

Theorem 6.3. Let m > 1. Suppose that v > 0 is a viscosity superso¬ 
lution in Qt with initial values v{x, 0) = 0 in If 

min{n-, j} G L^{0, T; 2 = 1,2,..., 

then 

(i) V G L'^(f2Ti) whenever q < 1^ and Ti < T, 

(ii) the function has a Sobolev gradient V(v^) G L'^'(Uti) when¬ 
ever q' < 1 + and Ti < T. 

The summability exponents are sharp. 

Proof. See [TT], Theorem 4.7 and 4.8. □ 

We start from the intrinsic Harnack inequality given in [3l Theorem 
3]. This is the fundamental analytic tool here. 

Lemma 6.4 (Harnack’s inequalty). Let m > 1. There are constants 
C and 7 , depending only on n and m, such that if u > 0 is a continues 
weak solution in 

CR^ 


i?(xo, 4i?) X (to — 40, to + 40), where 0 


n(xo,to)™“^ 
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then the inequality 


(6.5) 

is valid. 


u{xo,tQ)<^ inf u{x,tQ + 6) 

Br{xo) 


Again the waiting time 9 depends on the solution itself. Then we 
need the separable solution 

0(a;) 

{t - to)^ 

where the function S'” G Wl’‘^(yL) is a weak solution of the auxiliary 
equation 

A(0-) + —^ = 0, 
m — 1 

which is the Euler-Lagrange Equation of the variational integral 

/^|'u|"*+^dx 

This function is known as “the Friendly Giant”, see [V, p. Ill] and often 
serves as a minorant. When extended as 0 when t < to it becomes a 
viscosity supersolution in the whole hi x M. 


Proposition 6.6. Suppose that we have an increasing sequence 0 < 
hi < h 2 < ho <... of viscosity supersolutions in VLt and denote 
h = imik^oohk. If there is a sequence {xk,tk) —)■ (xo,to) such that 
hk{xk,tk) —)■ oo, where Xq E fl and 0 <to <T, then 

liminf h{y, t){t — to) ’"-i >0 for all x E fl. 

t>to 

Thus, at time to, 

lim h{y,t) = oo in hi. 

t>to 

Remark 6.7. Notice that the limit function is not a solution in the 
whole domain, but it may, nonetheless, be finite at each point. (This 
is different from the Heat Equation, see [23].) 

If it so happens that the subsequence in the Proposition does not 
exist, then we have the normal situation with a solution; 


Proposition 6.8. Suppose that we have an increasing sequence 0 < 
hi < h 2 < ho <... of viscosity supers olutions in Qt o,nd denote 
h = Ymik^oohk- If the sequence {hk} is locally hounded, then the limit 
function h is a supersolution in Qt- 

After this, the proof proceeds along the same lines as for the p- 
parabolic equation. A difference is that the infimal convolution should 
be replaced by the solution to an obstacle problem as in Chapter 5 of 

m 
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